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Abstract 

We study the distortion of electron neutrino energy spectrum due to oscillations 
with the sterile neutrino <-> i^s, for different initial populations of the sterile state 
5Ns at the onset of oscillations. The influence of this spectrum distortion on Big 
Bang Nucleosynthesis is analyzed. Only the case of an initially empty sterile state 
was studied in previous publications. 

The primordial abundance of He-4 is calculated for all possible 5Ng: < SNg < 1 
in the model of oscillations, effective after electron neutrino decoupling, for which 
the spectrum distortion effects on the neutron-proton transitions are the strongest. 

It is found that the spectrum distortion effect may be the dominant one not 
only in the case of small 5Ns, but also in the case of big initial population of Ug. 
For example, in the resonant case it may play a considerable role even for very big 
6N, ~ 0.8. 
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Introduction 



Sterile neutrinos Vg ™ay be present at the onset of the nucleosynthesis epoch. There 
may be different reasons for their production — they are naturally produced in GUT 
models [1] , in models with large extra dimensions [2] and Manyfold Universe models [3] , 
in mirror matter models [4]. They may be produced also in f^^,- <-> Vg oscillations in the 
preceeding epoch ^ according to the most popular 4- neutrino mixing schemes [5]. 

The degree of population oiusi ^^s, nia-y be different depending on the concrete model 
of i/g production. Hence, we will consider further on SNg as a parameter. 

In recent years, strong constraints on the sterile neutrino impact in oscillations, ex- 
plaining atmospheric and solar neutrino anomalies, were obtained from the analysis of 
experimental oscillations data [6, 7, 8]. And although a pure sterile channel solution is 
excluded for any of the neutrino anomalies, the analysis of neutrino oscillation experi- 
mental data allows a certain, small fraction of sterile neutrino to participate into these 
oscillations. In the analysis [9] and [10] it was suggested that a small fraction of Us is not 
only allowed but even desirable for solar neutrino data. Recent measurements of the cross 
section of '^Be{p,^)^B reaction lead to a predicted total neutrino flux by 13% larger 
than measured by SNO, which also may be providing evidence for sterile neutrino [11]. 

There also exist stringent cosmological constraints on i/g produced in oscillations, based 
on oscillations influence on BBN nucleosynthesis of "^He. For the case z/^ ^- ^ Ug sec 
for example refs. [12, 13] and for the z/g ^ Ug see refs. [14, 13]. The constraints on 
electron-sterile neutrino oscillations [15, 16], excluded the active-sterile LOW solution to 
the solar neutrino puzzle, in addition to the already excluded in pioneer works (sec for 
example ref.[17]) sterile LMA solution and sterile solution to the atmospheric neutrino 
anomaly. (For more detail about constraints on neutrino oscillations from cosmology, see 
refs. [18, 12].) 

Most of the cosmological constraints on active-sterile mixing were obtained in simple 
two-neutrino mixing schemes ^, for the case when the sterile neutrino state was initially 

^For example, atmospheric and LSND neutrino data require oscillations with maximal mixings and 
mass differences ^m„^„ ^ 10~^ eV^ and Smj^gj^j^ 10~^ eV^, which are effective before the BBN epoch. 
So, in many schemes owing to Ug or 1/^ ^ Vg oscillations, Vg state may be partially thermalized 

before the nucleosynthesis epoch. 

^For discussion and calculation of cosmological constraints in a specific 4-neutrino mixing schemes see 
ref. [19, 20, 13] 
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empty at the epoch before oscillations became effective in the Universe evolution, 6Ns = 
— 3 = 0. Ni, is the number of neutrino species in equilibrium. Since the presence 
of a non-empty sterile state before oscillations was not considered in previous analysis 
of oscillations effects on the neutrino spectrum distortion and on BBN, in this work we 
address this question. We omit the assumption SNs = and explore the general case 
SN, + 0. 

In the v^J^,^T <-> Vg oscillation case, bNg 7^ present initially just leads to an earlier 
increase of the total energy density of the Universe, and it is straightforward to re-scale 
the existing constraints. In the v^. Vg oscillations case, however, the presence of Vg at 
the onset of oscillations influences also the kinetic effects of i^e <-> Vg on BBN. Therefore, 
we chose <-> Vg oscillations case for exploring electron neutrino distortion, caused by 
oscillations, and its influence on nucleons freezing and on primordial ^le, Y.p for different 
bNg values in the range < 5A^s < 1. 

The spectrum distortion of the electron neutrinos 

We analyze the case of oscillations effective after neutrino decoupling, therefore further on 
we denote by ^Ng the degree of population of the sterile neutrino state at active neutrino 
decouphng (T ~ 2 MeV) . Here we explore f e spectrum distortion , considering the degree 
of population of the sterile neutrino state bNg as a free parameter and, varying its value 
in the range [0, 1] with a step 0.1. In the next section we calculate the kinetic effect of 
oscillations on primordial abundance of ^le for different ^Ng. 

The mixing of electron neutrino with Vg has the following two types of effects on BBN: 

(a) it leads to an increase of the energy density of the Universe [21], and 

(b) it changes the nucleons kinetics essential for n/j^-frcezing, through the depletion of 
electron neutrino number density, distortion of the equilibrium spectrum of Pg, and pro- 
duction of asymmetry between neutrinos and anti-neutrinos, all due to oscillations [18]. 
We will parametrize these kinetic effects and denote them SNkin further on. 

(a) The first effect is usually described by an increase of the effective number of the 
energy density degrees of freedom g^ff = (30/7r^)(p/T^). At the BBN epoch g^ff = 
10.75 -|- 7/46Ns{T,y/Ty. Hence, this effect leads to a faster expansion rate H ~ (jfgj^ and 
higher freezing temperature for nucleons Tf ~ 5*6// > when nucleons were more abundant: 

n/p ~ exp [—Am/Tf] 

This reflects into an overproduction of "lie, since it strongly depends on the n/p- 
freezing ratio: Yp ~ 2exp(— Am/T/)/[l -|- exp(— Am/T/)], where Am = rUn — rup ~ 1.3 
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MeV is the neutron-proton mass difference. This effect of gejf increase on BBN is well 
known [22]. The approximate fit to the exact calculations is: 6Yp ~ O.OlS^A^s. The 
maximum helium overproduction corresponding to 6Ns = f is ~ 5%. 

(b) The influence of the Icinetic effects of oscillations on BBN is quite obvious, hav- 
ing in mind that: (i) oscillations take place between equilibrium electron neutrino and 
less populated sterile neutrino ensemble, that (ii) the oscillations probability is inversely 
proportional to the energy of neutrinos P ~ 5rn? / so that neutrinos with different mo- 
menta start oscillating at different cosmic times, and that (iii) the proton density is bigger 
than the neutron one. Due to that, the neutrino energy spectrum ni,{E) may strongly 
deviate from its equilibrium form [23, 18]. 

In case oscillations proceed after the decoupling of active neutrinos, a strong spectrum 
distortion for both the electron neutrino and the anti-neutrino is possible. This spectrum 
distortion affects the kinetics of nucleons freezing - it leads to an earlier n/p-freezing and 
an overproduction of ^ie yield. 

The effect can be easily understood having in mind that the distortion leads both to 
a depletion of the active neutrino number densities in favor of the sterile ones N^: 



and to a decrease of the mean neutrino energy. ^ This lowers the weak rates, governing 
nucleons transitions during neutrons freezing, with respect to their values in the standard 
BBN model, F^^eait ~ N^^E^, and hence reflects into earlier freezing when neutrons were 
more abundant. So, helium is over-produced. 

The generation of neutrino-antineutrino asymmetry in the resonant oscillations has a 
subdominant effect on "lie. It slightly suppresses oscillations at small mass differences, 
leading to a decrease of helium overproduction. 

SNkin depend strongly on the initial population of the sterile neutrino at BBN. Larger 
6Ns decreases the kinetic effects, because the element of initial non-equilibrium between 
the active and the sterile states is less expressed. Hence, for any specific value of SNg it 
is necessary to provide a separate analysis. 

In the case SNg = 1 i/g are in equilibrium (the sterile state is as abundant as the 

electron one), and hence the n-p kinetics does not feel the oscillations, dNkin — 0. The 

final effect is only due to the energy increase, i.e. SNtot — 1- 

^The decrease of the electron neutrino energy due to oscillations into low temperature sterile neutrinos, 
has also an additional effect: Due to the threshold of the reaction converting protons into neutrons, when 
neutrinos have lower energy than the threshold one, protons are preferably produced, which may lead to 
an under-production of '^He [24]. However, this turns to be a minor effect in the discussed oscillations 
model. 



V 
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In the case 6Ns = the kinetic effect of oscillations was studied numerically for both 
the resonant [15] and non-resonant [16] oscillation cases. In this case SNkin for given 
fixed mixing parameters reach their highest value, 5A^^"^, as far as the non-equilibrium 
element — the difference between the sterile and active neutrino number densities at 
the beginning of oscillations — is the greatest. The overproduction of ^e may be enor- 
mous: up to 13.2% in the non-resonant oscillation case and up to 31.8% in the resonant 
one [25]. This corresponds effectively to a httle more than 6 additional neutrino states 
^^klT ~ 6. For the case SNg — the kinetic effects were carefully studied and strin- 
gent cosmological constraints on the oscillation parameters were obtained on the basis of 
helium-overproduction, caused by the kinetic effects of oscillations [14] . 

In this work, accounting for all the effects (a) and (b), we calculate Yp{5Ns, 5w?, sin^29) 
for < 5Ns < 1 values, and reveal the dependence of SN^in on SNg. 

We have analyzed the self-consistent evolution of the oscillating neutrino and the 
nucleons from the neutrino decouphng epoch at ~ 2 MeV till the freezing of nucleons. 
We have followed the line of work described in detail in ref. [15], omitting the assumption 
for negligible density of the sterile neutrinos at the onset of ^ oscillations. 

It is hardly possible to describe analytically, without some radical approximations, the 
non-equilibrium picture of active-sterile neutrino oscillations, producing non-equilibrium 
neutrino number densities and distorting the neutrino spectrum. Satisfactory precise 
analytical description was found only for the case of relatively fast oscillations proced- 
ing before neutrino freezing, with 5m? > 10~^ eV^ and small mixing angles [26, 13]. 
Therefore, we have provided a self-consistent numerical analysis of the evolution of the 
nucleons number densities n„ and the ones of the oscillating neutrinos p and p in the 
high-temperature Universe, using the following coupled integro-diffcrential equations, the 
first equation describing the kinetics of the neutrino ensembles in terms of the density 
matrix of neutrino p and anti-neutrino p, the second equation ~ the kinetic evolution of 
the neutrons. 



dpjt) dpjt) 



dt " dpu 

+i [Ho, p{t)] + iV2Gf {±C - Q/M^) [a, p{t)] + O (gI) , 



(1) 




(2) 
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where aij = U*^Uje, Pu is the momentum of electron neutrino, n stands for the number 
density of the interacting particles, dl](i, j, k) is a phase-space factor, and A is the ampli- 
tude of the corresponding process. The plus sign in front of C corresponds to the neutrino 
ensemble, the minus sign - to the anti-neutrino ensemble. 

Mixing just in the electron sector is assumed: i/j = Uu vi {I — e,s). The initial 
condition for the neutrino ensembles in the interaction basis is of the form: 

where n^'* = exp(— i?j^/T)/(l -|- exp{—E^/T)), while S measures the degree of population 
of the sterile state. 

Ho is the free neutrino Hamiltonian. The 'non-local' term Q arises as a W/Z propa- 
gator effect, Q El, T . C is proportional to the fermion asymmetry of the plasma and 
is essentially expressed through the neutrino asymmetries C ~ ^L^,^ -\- Lj,^ -\- L^,^, where 

L^,,r ~ {N^,r - N-^,f)/N^ and L,^ ~ ! ^Mpll - Pll)/N^- 

The equations are for the neutrino and neutron number densities in momentum space, 
which allows to describe precisely the kinetic effects: spectrum distortion and neutrino- 
antineutrino asymmetry growth due to oscillations. For the description of the spectrum 
1000 bins were used in the nonresonant oscillations case, and at least 5000 in the resonant 
one. These equations provide a simultaneous account of the different competing processes, 
namely: neutrino oscillations. Universe expansion, neutrino forward scattering, nucleons 
t r ansf ormat ions . 

The analysis was performed for all mixing angles d and mass differences (5m^ < 10"'^ 
eV^. The analyzed temperature interval was [2.0,0.3] MeV, because at temperatures 
higher than 2 MeV the deviations from the standard BBN model without oscillations are 
negligible in the discussed model of oscillations. 

As expected, the spectrum distortion is less expressed when increasing the degree of 
population of the sterile neutrino state 5Ns (Fig.l.). Correspondingly, the kinetic effect 
on primordial nucleosynthesis should decrease. The results of our numerical analysis on 
spectrum distortion at different 5Ns are illustrated in the Figs, la-c, where the dependence 
of the energy spectrum distortion of the electron neutrino on the initial population of the 
sterile state is shown. 

At each temperature we have plotted the spectrum for three different levels of initial 
population of the sterile neutrino, namely SNg = 0.0, 0.5, 0.8. 
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Figure la: The figure illustrates the degree of distortion of the electron neutrino energy spec- 
trum x^pLhix), where x = E/T at a characteristic temperature 1 McV, caused by resonant 
oscillations with a mass difference 5m? = 10~^ eV^ and sin^ 2'd = 0.1 for different initial sterile 
neutrino populations, correspondingly SNg = (the lower curve), SNg = 0.5 and dNg = 0.8 (the 
upper curve). The dashed curve gives the equilibrium neutrino spectrum for comparison. 
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Figures lb,c: Distortion of the electron neutrino energy spectrum at a temperature 0.7 MeV 
(Fig. lb) and 0.5 MeV (Fig.lc) for the same parameters as for Fig. la. 
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The oscillations parameters are \6m'^\ = 10~^ eV^ and sin^2'i? = 0.1. For illustrating 
the evolution of the spectrum distortion we have presented it at characteristic tempera- 
tures 1 (Fig. la), 0.7 (Fig.lb) and 0.5 MeV (Fig.lc). At each 6Ns the characteristic be- 
havior of the spectrum distortion due to oscillations is observed. Namely, since oscillation 
rate is energy dependent F ~ Sm^/ E the low energy part of the spectrum is distorted first 
(as far as low energy neutrinos start to oscillate first) and later the distortion penetrates 
noticeably into the more energetic part of the spectrum. 

We have found that the neutrino energy spectrum n^^E) may strongly deviate from its 
equihbrium form during all the period of interest (2 MeV - 0.3 MeV) even for considerably 
large 5Nsi and hence the spectrum distortion may constitute the dominant effect on the 
overproduction of ^ie. 

The kinetic effect 

For different SNg we calculate precisely the n/p- freezing, essential for the production of 
helium, down to temperature 0.3 MeV. Then we calculate Yp, accounting adiabatically 
for the following decay of neutrons till the start of nuclear reactions, at about 0.1 MeV. 

We have found that neutrino spectrum distortion effect on BBN is very strong even 
when there is a considerable population of the sterile neutrino state before the beginning 
of the electron-sterile oscillations. It always gives positive SNkin, which for a large range 
of initial sterile population values, are bigger than 1. The kinetic effects are the strongest 
for SNs = 0: Y^'''^{SN„Sm'^,sm^2i)) = ^^(0, ^m^, sin^ 2i?). They disappear for SN, = 1, 
when i/g and states are in equilibrium, and the total effect reduces to the SBBN with 
an additional neutrino. 

In Fig. 2 we present the frozen neutron number density relative to nucleons — 
Nl/Nnuc as a function of the sterile neutrino content at neutrino decoupling for a resonant 
and a nonresonant oscillation case. The oscillation parameters are Sm? = 10~^ eV^ and 
Sm^ = -10"^ eV^ and sin^ 2^ = 10"^ As far as SYp/Yp = SX^/X^, the results are 
representative for the overproduction of primordially produced helium. 
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Figure 2: The solid curves present frozen neutron number density relative to nucleons = 
N.[/Nnuc as a function of the sterile neutrino initial population. The dashed curves present only 
the kinetic effect, while the dotted curve presents the effect due to the energy density increase. 
The upper two curves (dashed and solid) correspond to the resonant case, the lower dashed and 
solid curves - to the nonresonant one. 



The dotted curve presents only the effect (a), due to the energy density increase 
Xl = f{SNs), the dashed curves present the pure kinetic effects (b) = f{5Nkin), while 
the solid lines give the total effect. The upper dashed and solid curves correspond to the 
resonant case, the lower ones to the non-resonant one. 

The analysis for these concrete oscillation parameters, shows that the overproduction 
of helium is strongly suppressed with the increase of SN^ for the resonant case, while in the 
non-resonant case it increases with SNg. This is a result of the fact that, in the resonant 
case, the kinetic effects (b) due to the spectrum distortion are the dominant contribution 
to the overproduction of helium, even for very large degree of population of the sterile 
state, while in the non-resonant case the main contribution comes from the increase of 
degrees of freedom already at very small ^A^^. An empirical approximation formula is: 

5Yp = 0.013[<5iV,^r(l - ^Ns) + 5Nsl 
where SNJ^^^ is the value calculated in the case of oscillations with an initially empty 
sterile state, i.e. SNf^t = ^^kin^i^ ~ ^^s) + ^^s- It is a good approximation for the 
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non-resonant case and a rather rough one for the resonant case: the deviation from the 
exactly calculated helium given in Fig. 2 may be up to 5Yp/Yp ~ 0.8%. Still, it can give 
some idea of 5Yp/Yp dependence on 5Ns. 

However, for other mixing parameters, the kinetic oscillation effects in the non-resonant 
case can be also considerable, as shown in ref.[25], the kinetic effect can be as high as 
^ 3 for initially empty sterile state. Hence, in the non-resonant case the spectrum 
distortion effects may be the dominant one even for much larger 8Ns than in the case 
illustrated in Fig. 2. 

For each concrete SNg value a detailed numerical analysis is necessary to reveal the 
interplay of effects (a) and (b) and their influence on primordial production of ^ie. 

In a forthcoming paper we apply the results obtained here to define the isohelium 
contours corresponding to 3% overproduction of ^ie for different SNg, and we present the 
cosmological constraints for nonzero 5Ns. 

Conclusions 

The presence of a non-empty sterile state before <-> oscillations was not considered 
in previous analysis of oscillation effects on the neutrino spectrum distortion and 

on BBN. In this work we have studied the kinetic effects due to <-> Ug oscillations in 
the general case Q < SNg <1. 

We have provided a numerical analysis, investigating how the presence of the sterile 
neutrino state, partially populated before oscillations, will influence the production of 
^le in the model of BBN with electron-sterile oscillations effective after electron neutrino 
decoupling. 

We have found that the effect of the neutrino spectrum distortion due to oscillations 
may he very strong, even for a considerable initial population of the sterile neutrino state. 
Correspondingly, the kinetic effect of oscillations remain the dominant one even for big 

The results of this analysis may be applied for different models generating sterile 
neutrino, like GUT models, mirror models, extra-dimensions models, etc., as far as the 
initial value of population of the sterile state 5Ns depends on the concrete model of its 
production. These results may be of interest also for mixing schemes in which a portion 
of Vg have been brought into equilibrium before neutrino decoupling, due to z/^ ^ Vs 
or Vr ^ oscillations. In case the Vs presence is due to the much earlier (at atmo- 
spheric mass difference scale, or LSND) oscillations of z/^ ^- ^ z/^, 5Ns may be directly 
connected with the available constraints on the sterile neutrino fraction, deduced from 
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the neutrino oscillations experimental data analysis. So, we hope that the results may be 
indicative and helpful for choosing among the different possibilities for the sterile fraction 
in the subdominant active-sterile oscillations used in the oscillation analysis of neutrino 
anomalies. 

A more general study of kinetic oscillations effects on BBN for non-empty initially 
sterile state in the framework of 4- neutrino mixing schemes seems appropriate, although 
much more complicated. 

I thank A. Dolgov for useful comments on the draft version, M. Chizhov and A. 
Strumia for stimulating discussions during the preparation of this work. I am grateful to 
my parents for their help with the kids, during the work on this paper. 

I appreciate the Corresponding Associate position at CERN TH and the Regular 
Associateship at the Abdus Salam ICTP, Trieste. This work was supported in part by 
Belgian Federal Government (office for scientific affairs grant and lAP 5/27). 

References 

[1] Brahmachari B., Phys. Lett. B461, 243 (1999); Europhys. Lett. 46, 428 (1999); Mo- 
hapatra R., in "Current Aspects of Neutrino Physics" ed. D. Caldwell (Springer- 
Verlag), hep-ph/9910365 (1999); B. Brahmachari, and R. Mohapatra, Phys. Lett. 
B437, 100 (1998); Mohapatra R., Nandi S. and Perez-Lorenzana A., Phys. Lett. 
B466, 115 (1999); Mohapatra R., Nucl. Phys. Proc. Suppl. 91, 313 (2000). 

[2] Dvah G., A. Smirnov, Nucl. Phys. B563, 63 (1999); Mohapatra R., S. Nandi, A. 
Prez-Lorenzana, Phys. Lett. B466, 115 (1999); loannisian A., Pilaftsis A., Int. J. 
Mod. Phys. A 16, 931 (2001); loannisian A., and J. W. F. Valle, Phys. Rev. D 63, 
073002 (2001). 

[3] Arkani-Hamed N. et al, JHEP 0012, 010 (2000). 

[4] Li T., Yang C., Phys. Rev. 104, 254 (1956); Kobzarev Y., L Okun, I. Pomeranchuk, 
Yad. Fiz. 3, 1154 (1966); Blinnikov S., Khlopov M., Yad. Fiz. 36, 675 (1982); Foot 
R., H. Lew, R. Volkas, Phys. Lett. B 272, 67 (1991); Ahmedov E., Berezhiani Z., 
Senjanovic G., Phys. Rev. Lett. 69, 3013 (1992); Berezhiani Z., Mohapatra R., Phys. 
Rev. D52, 6607 (1995); Berezhiani Z., Dolgov A., Mohapatra, R., Phys. Lett. B 375, 
26 (1996); Berezinsky V., Narayan M., Vissani F., Nucl. Phys. B658, 254 (2003). 



12 



[5] Barger V., T. J. Weilcr, K. Whisnant, Phys. Lett. B427, 97 (1998); Bilenky S., G. 
Giunti, W. Grimus, Nucl. Phys. S. B 77, 151 (1999); Valle J., Phys. Atom. Nucl. 
63 (2000); Barger V. ct al., Phys. Lett. B489, 345 (2000); Mohapatra R., in Current 
aspects of neutrino physics, ed. Caldwell D., p. 217; Guinti C. and Laveder M., JHEP 
0102, 1 (2001); Fogli G., E. Lisi, A. Marrone, Phys. Rev. D 64, 093005 (2001); 
Maltoni M., Schwetz T., Valle J., Phys. Lett. B518, 252 (2001); Peres O., Smirnov 
A., Nucl. Phys. B599, 3 (2001); Maltony M., Schvetz T., Valle J., Phys. Rev. D 65, 
093004 (2002); Bahcall J., Pena-Garay C., JHEP 0311, 004 (2003); Maltoni M. et 
al., hep-ph/0305312. 

[6] Gonzalez- Garcia M., and Pena-Garay C., Nucl. Phys. B91, 80 (2001); Grimus W., 
Acta Phys. Polon. B30, 3067 (1999); Maltoni M. et al, Nucl. Phys. B643, 321 (2002); 
Barger V. et al., Phys. Lett. B537, 179 (2002); Barger V., D. Marfatia, K. Whisnant, 
Phys. Rev. Lett. 88, 011302 (2002); Giunti C., 2003, hep-ph/0309024; Maltoni M. et 
al, Phys. Rev. D 67, 013011 (2003); Maltony M., Schwetz T., Valle J.: 2003, Phys. 
Rev. D 67, 093003 (2003). 

[7] Strumia A., Phys. Lett. B539, 91 (2002); Nakaya T., hep-ex/0209036, 2002. 

[8] Bahcall J., Gonzalez-Garcia M.C., and Pena-Garay C., JHEP 0207, 054 (2002); 
Bahcall J., Gonzalez-Garcia M., and Pena-Garay C., JHEP 0302, 009 (2003); Bahcall 
J., Pena-Garay C., JHEP 0311, 004 (2003). 

[9] Barger V. et al, Phys. Lett. B 489, 345 (2000). 

[10] Holanda P., Smirnov A., 2003, hep-ph/0307266. 

[11] Gai M., nucl-ex/0312003, 2003. 

[12] Dolgov A., Phys. Rept. 370, 333 (2002). 

[13] Dolgov A., Villante F., 2003, hep-ph/0308083. 

[14] Kirilova D. P., and Chizhov M. V., Nucl. Phys. S. B 100 360 (2001). 

[15] Kirilova D. P., and Chizhov M. V., Nucl. Phys. B 591, 457 (2000). 

[16] Kirilova D. P., and Chizhov M. V., Phys. Rev. D 58, 073004 (1998). 

[17] Barbieri R., and Dolgov A., Phys. Lett. B 237, 440 (1990); Nucl. Phys. B 349, 743 
(1991). 



13 



[18] Kirilova D. P., and Chizhov M. V., Phys. Lett. B 393, 375 (1997); CERN-TH/2001- 
020, astro-ph/0108341. 

[19] Bilenky S. et al., Astroparticle Phys. 11, 413 (1999); Okada N., and Yasuda O., Int. 
J. Mod. Phys A 12, 3669 (1997). 

[20] Di Bari P., Phys. Rev. D65, 043509 (2002); Di Bari P., astro-ph/0302433. 

[21] Dolgov A., Sov. J. Nucl. Phys. 33, 700 (1981). 

[22] Hoyle F., and Tayler R., Nature 203, 1108 (1964); Shvartsman V., JETP Lett. 9, 
184 (1969). 

[23] Kirilova, D. P., 1988, preprint JINR E2-88-301. 

[24] Dolgov A.,Kirilova D., Int. J. Mod. Phys. A3, 267 (1988). 

[25] Kirilova D. P., Astropart. Phys 19, 409 (2003). 

[26] Dolgov A. D., 2003, liep-ph/0306154. 

[27] Dolgov A. et al., Nucl. Phys B632, 363 (2002). 

[28] Trotta R., Hansen S., 2003, astro-ph/0306588. 



14 



